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\ Abstract. A new class of nonlinear partial differential equations with distributed 

\ in space and time state-dependent delay is investigated. We find appropriate assump- 
tions on the kernel function which represents the state-dependent delay and discuss 
advantages of this class. Local and long-time asymptotic properties, including the 
existence of global attractor, are studied. 
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1. Introduction 

Theory of delay differential equations is one of the oldest and simultaneously, in- 
tensively developing branches of the theory of infinite-dimensional dynamical systems. 
This theory covers ordinary and partial delay differential equations, includes studies 
of discrete and distributed, finite and infinite delays. Classical methods of differen- 
tial equations, theory of distributions and functional analysis allow one to study wide 
classes of ordinary and partial differential equations with delay. We mention only sev- 
eral monographs which are classical references for delay equations [91 [7J [lj [HI [31] and 
also works which are close to this investigation [281 El El EEE1 El [201 EI] • Nevertheless, 
each (nonlinear) equation requires a separate and careful studying. 

Recently, a new class of delay equations attracts attention of many researchers. 
These equations have delay (delay term) which may change, according to the state of 
the system i.e. state-dependent (state-selective) delay. The study of such equations 
was started in the case of ordinary equations [151 [HI EQ] and it was recently continued 
for P.D.E.s in [20]. [21] . For more detailed discussion and references on delay equations 
see e.g. introduction in [2D]. We continue our previous research J2DJ, [2T] and present a 
wider class of nonlinear equations with distributed in space and time state-dependent 
delay terms. 

Let us illustrate the main question studied in this article on the simplified object 
which is a local in space delay term. Consider the following simple distributed in 
time delay term J® r b(u(t + 6,x))£(0)d6. Here function £ belongs to some space of 
real valued functions defined on the delay interval (— r, 0). This (kernel) function £ 
represents the rule how the information on the previous stages of the system (function 
u) is used to model the process. As discussed (see e.g. [201 [21]), this rule may change 
according to the state of the system. Let us denote by v G H the state coordinate, 
where H represents the phase space. With this notations the state-dependent (state- 
selective) delay rule reads v) : (— r, 0) x H — ► R and the corresponding delay term 
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becomes J_ r b(u(t + 9,x))£(9,v)d9. This is a simplified (local) example of delay terms 
studied in [201 [21]. As we will see (section 3), studying some questions (e.g. stationary 
solutions), there is a need to use a wider class of functions £ (delay rules) which are 
space-dependent i.e. £(9,x,v). For example, considering a biological system, where 
u(t, x) represents the density of a population at time moment t at point xGll, the 
delay rule £(0,v) is the same for all points x in the domain Q, while the delay rule 
£(9, x, v) depends on the point x G Q (e.g. due to the dependence of food resources on 
points in fl). This interpretation shows that the delay rule £(9,x, v) is more realistic 
biologically and as we will see (section 3) mathematically. Taking into account the 
above motivation we need to find an appropriate class of functions £, concentrating on 
the character of dependence of £ on the coordinate igSI. This is the main goal of the 
article. It is interesting to mention that in spite of the fact that for any time moment 
t > solutions u(t) belong to the space L 2 (Q) and the phase coordinate (u(t); u(t+9)) 
belongs to L 2 (Q) x L 2 (— r, 0; L 2 (Q)), the values of function £ (as functions of x) do 
not necessary belong to L 2 {Vl). They belong to a wider space D(A^ 1 ^ 2 ) D L 2 {Vl) (for 
more details see theorems 1 and 2 below). 

The proposed model has an essential advantage in comparison with the previous 
ones (see [201 EI]) to cover the case of finite and even infinite sequences of isolated 
stationary solutions. We also present an algorithm to construct such state-dependent 
delay terms. 

The article is organized as follows. In section 2 we present the model, prove the 
existence and uniqueness of weak solutions, construct the dynamical system and prove 
the existence of a global attractor. Section 3 is devoted to stationary solutions and 
the possibility to use our system to construct a dynamical system with an a-priory 
given set of isolated stationary solutions. The results may be applied to the diffusive 
Nicholson's blowflies equation. 

2. Formulation of the model with distributed delay 

Consider the following non-local partial differential equation with state- dependent 
distributed in space and time delay 

■§tu(t, x) + Au(t, x) + du{t, x) 
= I°r {In Ku(t + 9, y))f(x - y)dy} £(9, x, u(t),u t )d9 (1) 
= (F(u t ))(x),x G n, 

where A is a densely-defined self-adjoint positive linear operator with domain D(A) C 
L 2 (Q) and with compact resolvent, so A : D(A) — ► L 2 (Q) generates an analytic 
semigroup, Q is a smooth bounded domain in R n °, f : Q — Q ^ R is a bounded 
function to be specified later, b : R — > R is a locally Lipschitz bounded map < 
Cb with Cb > 0), d is a positive constant. As usually for delay systems (see [9]) for 
any function u(t),t G [a, b], b > a + r with values in a Banach space X, we denote by 
Ut = u t (9) = u(t + 9), which is a function of 9 G [— r, 0] with parameter t G [a + r, b\. 
Constant r > is the (maximal) delay of the system. 

The function £(•, -, •) : [— r, 0] x Q x H — >• R represents the state-dependent dis- 
tributed delay. We denote for short H = L 2 {VL) x L 2 {— r, 0; L 2 (Q)) and also use || ■ || 
and (-, •) to denote the norm and scalar product in L 2 (Q). 
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We consider equations ([T]) with the following initial conditions 

u(0+) = u° E L 2 (n), u| ( _ r)0) = <p G L 2 (-r,0;L 2 (ft)). (2) 

So we write (n°, ip) G if. 

Now we study the existence and properties of solutions for distributed delay prob- 
lem CO), ©. 

Definition 1. A function u is a weak solution of problem (TJP subject to the ini- 
tial conditions (TJj on an interval [0,T] if u G L°°(0, T; L 2 (fi)) fl L 2 (—r, T; L 2 (Q)) n 
L 2 (0, T; u(0) = <p(6) for 9 G (-r, 0) and 

- [ T (u,v)dt+ [ T (A^u,Ah)dt+ [ T (du - F(u t ),v)dt = -(u°,v(0)) (3) 
Jo Jo Jo 

for any function v G L 2 (0, T; D(A?)) with v G L 2 (0, T; and u(T) = 0. 



Theorem 1. Assume that 

(i) b : R — > i? is locally Lipschitz and bounded i.e., there exists a constant Cb so 
that that \b(w)\ < for all w G R; 

(ii) f : IT^Tl ^Ris bounded (|/(-)| < M f ); 

(Hi) £ : [— r, 0] x Q x L 2 (Q) x L 2 (—r, 0; L 2 (fi)) — ► i? satisfies the following conditions: 
a) /or any M > taere exists so i/iaf /or all (v\ ip 1 ) G satisfying 

lbi| 2 + /-rll^(s)|| 2 rfs<M 2 ,i = l,2 oneaas 



|v 1 -t; 2 ll a + 



||^(s)-^ 2 (s)|| 2 ds 



d6 

,1/2 



b) t/jere exists Cr^im > so inai 

llf (0, • , «, IU(A-V2) d# < C( € _i/2) /or a// (u, ^) G F. 



(4) 



(5) 



Taen for any (n ,<£>) G if = L 2 (tt) x L 2 (— r, 0; L 2 (f2)) the problem (TJP subject to 
the initial conditions (TJ|) aas a u>ea£; solution u(t) on every given time interval [0, T] 
and tazs solution satisfies 

u(t)EC([0,T];L 2 (n)). (6) 



Remark. Properties (Hi) -a) and (iii)-b) mean that£ as a function of the third and 
fourth coordinate (v,if)) G H is a (nonlinear) locally Lipschitz and globally bounded 
mapping £ : H — > r, 0; _D(v4~2)). 
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Proof of Theorem 1. Let us denote by {e k }kLi an orthonormal basis of L 2 (Q) such 
that Ae k = X k e k , < Ai < . . . < A^ — ► +oo. We say that function u m (t,x) = 

m 

9k m{t)ek{x) is a Galerkin approximate solution of order m for the problem (QJ), p|) 

k=l 

if 

j (u m + Au m + du m - F(uf), e k ) = 0, 

\ (u m (0+), e k ) = (u°, e k ), (u m (6), e k ) = (<p(6), e k ), V0 G (-r, 0) 

Vk = 1, . . . , m. Here g k , m G C x (0, T] R) n £ 2 (-r, T; i2) with g k , m (t) being absolutely 
continuous. 

Equations (|7j) for fixed m can be rewritten as a system for the Tri- 
dimensional vector-function v(t) = v m (t) = (gi tm (t), ■ ■ ■ , g m , m (t)) T . We notice that 

m 

h m (t,-)\\h(fl)= Eslm® = WW**- 

The standard technique (see e.g. [9]) gives that for any initial data ip G 
L?{— r, 0; R m ), a G R m there exist a > and a unique solution of J7} w 6 
L 2 (— r, a; _R m ) such that f o = V 9 an d f(0) = a, and f|[o, Q ] G C([0, a]; R m ) (for more 
details see Theorem 6 and Remark 9 from |19| and also Lemma from I21H . 

It is easy to get from (jHJ) and the boundedness of b and / that 



\(F(u t ),v) L 2 (n) 

f-0 



b(u{t + 6,y))f{x-y)dy 



£(#, x, u(t),u t )d8 > v(x)dx 



b(u(t + 6,y))\ / f(x -y)£{0,x,u(t),u t )v(x)dx \dy 



dd 



< c b M f \n\ / u(e, ;u(t),u t )\\ DiA - 1/2) de ■ \\v\\ D{A1/2) . 

J—r 

Using (jSJ), one has 

\(F(u t ),v) L2(n) \ < C b M f \Q\C { ^ 1/2) ■ \\A^ 2 v\\. 



(8) 



Now, we will get an a-priori estimate for the Galerkin approximate solutions for 
the problem ([I]),©. We multiply ([7j) by g k m and sum over k — 1, • • • , m. Hence for 
u(t) = u m (t) and t G (0, a] = (0, a(m)], the local existence interval for u m (t), we get 



1 d 



2dt 

Using dSJ, © we obtain 



u(t)\\ 2 + \\A^u(t)\\ 2 + d\\u(t)\\ 2 <\(F(u t ),u(t))l 



d „ 

~r \\ u 
dt u 

Since |||«(t)|| 2 



2 + p 1/2 u(t)|| 2 + 2d\\u(t)\\ 2 < C 2 M 2 \Q\ 2 Cl_ 1/2) = ~h- 



^MtW + 2d||«(t)|p 



at 



(i 



u 



(9) 



(10) 



2 dr 



+2df t \\u(T)\\ 2 dr), we denote by X (t) = \\u(t)\\ 2 + J^ \\A 1 / 2 u(r)\\ 2 dr + 2d \\u(T)\\ 2 dT 
and rewrite the last estimate as follows 4fX{t) < AVe obtain \if) < x(0) + kit = 
||u(0)|| 2 + k\t. So, we have the a -priori estimate 



u{t)\\ 2 + / \\A 1/2 u(T)\\ 2 dr + 2d / ||M(r)|| 2 dr < ||u(0)|| 2 + ^t 



(11) 



Estimate (jllj) gives that, for u° G L 2 (Q) the family of approximate solu- 
tions {w m (t)}^ =1 is uniformly (with respect to m G N) bounded in the space 
L°°(0,T;L 2 (ft)) n L 2 (0,T;D(A 1 / 2 )), where D(A 1/2 ) is the domain of the operator 
A 1 / 2 and [0, T] is the local existence interval. From (ITT]) we also get the continuation 
of u m (t) on any interval, so ( fTTI) holds for all t > 0. 

Using the definition of Galerkin approximate solutions (J7|) and their property 
(fTTI) . we can integrate over [0,T] to obtain J T ||A- 1 / 2 M m (r)|| 2 rfr < C r for any T. 
These properties of the family {u m (t)}^ =1 give that {(u m (t); u m (t))}™ =1 is a bounded 
sequence in the space 

X T = L°°(0, T; L 2 (fi)) n L 2 (0, T; D(A 1/2 )) x L 2 (0, T; D(yT 1/2 )). (12) 

Then there exist a function (tt(t); w(t) ) and a subsequence {u mfe } C {u m } such that 

(u mk ;u mk ) *-weakly converges to (u;u) in the space X?. (13) 

By a standard argument (using the strong convergence u mk — > u in the space 
L 2 (0,T; L 2 (f2)) which follows from ( fl3l) and the Doubinskii's theorem, one can show 
(see e.g. Lions (1969), Chueshov (1999) and Rezounenko (1997)) that any *-weak 
limit is a solution of ([1]) subject to the initial conditions ([2]). To prove the continuity 
of weak solutions we use the well-known (see also [TBT thm. 1.3.1]) 

Proposition 1 (Proposition 1.2 in [23]). Let the Banach space V be dense and 
continuously embedded in the Hilbert space X; identify X = X* so that V •—>■ X 
V*. Then the Banach space W p (0,T) = {u G L p (0,T;V) : u G L q (0,T;V*)} (here 
p~ x + q~ x = 1) is contained in C([0, T]; X). 

In our case X = L 2 (fi), V = D(A 1 ^ 2 ), V* = D(A~ 1/2 ),p = q = 1/2 (see (UJ),^). 
Hence Proposition 1 gives ([6]). The proof of Theorem 1 is complete. ■ 

Now we describe a sufficient condition for the uniqueness of weak solutions. 

Theorem 2. Assume that functions b and f are as in Theorem 1 (satisfy properties 
(i),(ii)) ; function £ satisfies property (Hi) -a) and 

£(•, •, v, G L°°(-r, 0; D(A^ 1 ^ 2 )) for all (v, i/>) G H. (14) 

Then solution of (T7J) ; given by Theorem 1 is unique. 

Proof of Theorem 2. Let u 1 and u 2 be two solutions of ([TJ, ([2]). Below we denote 
for short w(t) = w m (t) = u 1,m (t) — u 2,m {t) - the difference of corresponding Galerkin 
approximate solutions. Hence 

^ll^(t)|| 2 + 2||^ 1 / 2 ^(t)|| 2 + 2rf||^(t)|| 2 = <^(^) - ^(^),^(t)>. (15) 
Let us consider the difference (F{u\) — F{u 2 ),w{t)) in details (see (OQ)). 

(F(ul)-F(u 2 t )Mt))= [ 



b{u\t + 9, y))f(x - y)dy £(0, x, u\t), u\)d9- 
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b(u 2 (t + 0, y))f{x - y)dy £(0, x, u 2 {t),u 2 t )d9 



w(t, x)dx 



r un 



b(u\t + 9, y))f(x - y)dy \ £(0, x, u\t),v^)d9 



r Un 



b{u\t + 0, y))f(x - y)dy £(0, x, u\t),u\)dO 



■ w(t, x)dx, 



b(u 2 (t + 0, y))f(x - y)dy £(0, x, u L (t),u l t )d6 



n 



-r Uf! 

b{u\t + 9, y))f(x - y)dy} £(0, x, u 2 (t),u 2 )d9 
Using the local Lipschitz property of b, flHJ) and (j3J, we deduce 

-r {X 



\{F{u])-F{u 2 )Mt))\<L h M f 



(t + 9,y)\dy jf |£(0, x, )| • |w(t, z) | dx} 

H-cyw/ini 1 1 e , • , ^ x (* ) , ^) - € , - , « 2 (*) , ) 1 1 ^ ( ) ^ - i I -4 1/2 ^ (*) i I 



de 



+C b M f \fl\L 



w{t)\\ 2 + / \\w(t + s)\\ 2 ds 



1/2 



< ^My^fij ess sup \\Z(e,;u l (t),u l t )\\ D(A _ 1/2] 



6»e(-r,0) 



P 1/2 ^WII 

w (t + 0,-)ll-p 1/2 w (t)llrf0 



< -\\A 1 / 2 w(t)\\ 2 +-L 2 b M 2 \Q\ 



\w(t)\\ 2 + I \\w(t + s)\\ 2 ds 

l2 



ess sup ||£(0, -,u (t),u t )\\ D(A _ 1/2) 

Oe(-r,0) y ' 



\w(t+9,-)\\ 2 d9 



+ l\\A 1/ MtW + lc 2 b M 2 M 2 Ll 



\w{t)\\< 



\w(t + s)\\ 2 ds 



Finally, we get the existence of positive constants C\,C 2 such that 

r-0 



dt 



| (F(ul) - F(u 2 ),w(t))\ < \\A^w{t) H 2 + Ci / |K* + 0) \\ 2 d9 + C 2 \\w(t) \\ 2 
The last estimate and ffT5l) give 

\w(t)\\ 2 +2\\A^ 2 w(t)\\ 2 +2d\\w(t)\\ 2 < \\A l / 2 w(t)\\ 2 +Ci f° \\w(t+9)\\ 2 d9+C 2 \\w(t)\\ 



KWA^wWlf + dl I \\w(9)\\ 2 d9+ [ \\w(s)\\ 2 ds) + C 2 \\w(t)\\ 2 . 



Hence 

d 
dt 



\w 



(t)\\ 2 +\\A 1/2 w(t)\\ 2 +2dMt)\\ 2 <d 



w(9)\\ 2 d9+ / \\w 



,(s)\\ 2 ds)+C 2 \\w(t)\f 



and property j | ^4. 1 / 2 -^ 1 1 2 > Ai||u || 2 gives 



d_ 

If 



w(t)\\ 2 + (\ 1 + 2d) j\\w(s)\\ 2 ds < d \\w(6)\\ 2 d6 + J* \\w(s)\\ 2 ds^+C 2 \\w 



It implies that there exists C 3 > 0, such that for Z(t) = (lu'(t)ll 2 + (Ai + 
2d) Jq ||w(s)|| 2 <is, we have 



d 
~dt 



l dQ. 



-Z(t)<C z Z(t) + C x J r \\w(6)\\< 
Gronwall lemma implies 

Z(t) < ('||w(0)|| 2 + C 1 C 3 - 1 j° \\w(9)\\ 2 de) -e Cst . (16) 
The last estimate allows one to apply the well-known 

Proposition 2. [321 Theorem 9] Let X be a Banach space. Then any *-weak 
convergent sequence {wk}^^ G X* *-weak converges to an element Woo G X* and 
\\woo\\x < liminf^oo \\w n \\ X - 

Hence, for the difference w 1 (t) — u 2 (t) of two solutions we have 

\\u\t) - u 2 (t)\\ 2 + 2(Ai + d) C \u\s) - u 2 (s)\\ 2 ds 

Jo 

< (\\u\0) - u 2 (0)\\ 2 + C1C3- 1 fjW\e) - V 2 {d)\?de) ■ e c *. (17) 

We notice that by ([6]) the difference ||u 1 (t) — w 2 (£)|| makes sense for all t G 
[0,T], VT > 0. The last estimate gives the uniqueness of solutions and completes 
the proof of Theorem 2. ■ 

Theorems 1 and 2 allow us to define the evolution semigroup St : H —>■ H , with 
H = L 2 {Sl) x L 2 (-r, 0; L 2 (fi)), by the formula S t (u°; <p) = (u(t); u(t+6)), 9 G (-r, 0), 
where u(t) is the weak solution of (JTJ) , (J2J) . The continuity of the semigroup with 
respect to time follows from ([6]), and with respect to initial conditions from (ITTj) . 

For the study of long-time asymptotic properties of the above evolution semigroup 
we recall (see e.g. [2| 127]) 

Definition 2. A global attractor of the semigroup S t is a closed bounded setU in H, 
strictly invariant (S t U = U for any t > 0), such that for any bounded set B C H we 
have lim sup{dist H (S t y, U), y G B} = 0. 

Theorem 3. Assume functions b and f satisfy properties (i), (ii) of Theorem 1. Let 
function £ satisfy properties (iii)-a) of Theorem 1, ( fi^| ) and also there exists C^q > 
such that (c.f. fifty) 

f° U(0, -,v,ip) II d6 < C m for all (v, i/>) G H. (18) 

J — r 
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Then the dynamical system (S t ; H) has a compact global attractor U which is 
a bounded set in the space H 1 = D(A a ) x W, where W — {<p : <p G 
L°°(-r,0;D(A a )),ipe L°°(— r, 0; D(A a ~ 1 ))}, a < \. 

Proof of Theorem 3. To prove the existence of the global attractor we use classical 
theorem saying that it is sufficient for the dynamical system (St, H) to be dissipative 
and asymptotically compact (see [2l [27| HJ). 

The property (|T8|) gives the estimate stronger than ([8]): 

\(F(u t ),v) L 2 {n) \ < C b M f \Q\C m ■ \\v\\ (19) 

which is necessary for the property of dissipativeness of (St; H). The rest of the proof, 
including the property of asymptotic compactness, is standard (see e.g. [2JIUEI8] an d 
also [201 HQ) • ■ 

3. Stationary solutions 

For simplicity of presentation, in this section we consider operator A = (—Ad) > 
0, where is the Laplace operator in L 2 (Q) with the Dirichlet boundary conditions. 
In this case (which is sufficient for the application to the Nicholson's blowfly equation), 
we have D(A) = H 2 (n) n H%(Q), D(A 1 ' 2 ) = Hq(VL), D(A^ 1 ^ 2 ) = H~\n). For more 
details on this classical Sobolev spaces see e.g. [T3] . 

In this section we concentrate on the stationary solutions. First of all, by defini- 
tion 1 (of a weak solution), u(t,x) G L 2 (0, T; D(A 1 ^ 2 )) = L 2 (0, T; Hq(Q)), so for the 
stationary solution u(t,x) = u st (x), one has u st G Hq(Q). 

Let us consider an arbitrary function u st G Hq(Q) C L 2 (Q). Our goal is to find 
conditions on a function £(•,*,-,*) such that the system has stationary solution 
u(t) = u st G H£(Q) for all t G R. Let us denote by = v*(0) = u st , 9 G [-r, 0]. 

Since for u st = G Hq(Q) we can choose £(-,-, 0, 0) = 0, we concentrate below on 
the case u st ^ G H^(Q). 

From ([!]) and ^u(t,x) = 0, we have 



Au s \x) + d ■ u s \x) = J n b{u st {y))f(x - y)dy ■ f° r £(6, x, u st , u st )d9, xeQ. 

(20) 

As we will show, it is sufficient to define in a proper way the value of £ for the second 
and third coordinates equal (u st ,u st ) G H = L 2 (fi) x L 2 (— r, 0; L 2 (fi)) only. We 
propose to look for this value, decomposing it on the time and space coordinates i.e. 

((M 1 « s ',« s )=x(fl)^W, 0e[-r,O], iGft. (21) 

Now equation fl2"Ul) reads 

AM st (a;) + d ■ w st (a;) = / n 6(u s *(y))/(2; - y)dy ■ v(x) ■ J° r x(d)d9, x G O. 

(22) 

We need the following elementary 

Lemma. Assume n st / 6 L 2 (fi). Let the function f be strictly positive and 
f G C°°(Q — Q). Let function b be bounded and satisfy b(w) > for all w ^ 0. 



S 



Then the function 



p(x) = / b(u st (y))f(x-y)dy 



(23) 



satisfies properties: p G C(f2), mf{p(x) : x G f2} = p m j„ > and sup{^-p(x) : x G 
H,i= l,...,n } ^K™* < oo. 

The continuity of p on follows immediately from the continuity of / and the 
Cauchy-Schwartz inequality 



^(a; 1 ) -p(x 2 )\ 



< \ \b(u 



St I 



b(u st (y)) [fix 1 -y)- f(x 2 - y)] dy 
f(x - y) - f{x - y) dy 



(24) 
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We also use that for all y G Vt one has ((a; 1 — y) — (x 2 — y)\ 

Properties b(w) > for all w^O, u st ^ G L 2 (f2) and strict positivity of / imply 
that p(x) > for all x G £1. Hence, the continuity of p and the Weierstrass theorem 
give ini{p(x) : x G Q} = p m in > 0. The boundedness of partial derivatives of p is due 
to / € C°°(fW2). I 

We also assume 



x(0)</0 ^ 0. 



(25) 



Under the assumptions of Lemma and (125]) we have p(x) ■ f° r x(0)d9 ^ for all 
iGfi. So we can write (see ([22])) 



Au (x) + d • u (x 



v\x) 



InKu st {y))f{x-y)dyfl r x{o)de 



(26) 



We notice that (126}) is the equality in .D(A -1 / 2 ) = H~ l (Vt) (in the sense of distribu- 
tions). As we saw, by definition 1 (of a weak solution), u(t, x) G L 2 (0,T; D(A l l 2 )) = 
L 2 (0,T; iifg (Q)) , so for the stationary solution u(t,x) = u st (x), one has u st G Hq(£1). 
This implies Au st G = D(yl _1 / 2 ) and Au st + d ■ u st E i? _1 (0). To show that 

t> G i/ _1 (n) we remind the following 

Proposition 3. [T3l Theorem 12.1] Let m be positive integer. Then any element 
h G H~ m (Q) may be represented (in the non-unique way) in the form 

h = D j hj, hj G L 2 {Q). 

\j\<m 



Here D c 



d a i- 



a 



{ai, . . . , a n }, 



a 



oti + . . . + a n . 



dx^...dxZ r ' 

In our case m = 1 and if we denote by h = Au st + d ■ u st G H~ l {Vt) and by 
q(x) = p _1 (x), then v G H~ l {VL) reads as qh G H^iVt). Using proposition 3, we 
write h = h + E?=i hi G L 2 {Vt). Hence, 



"0 Q n 

qh = qh + Y] q—K = qh + V 
r~f OXi f-f 



-^-(qhi) - K-^-q 
oxi axi 



(27) 
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Remark. We notice that all the derivatives are understood in the sense of 
distributions (see Uffy, \ffl$). The term q ■ h is understood as the distribution which is 
obtained by multiplication of the distribution h by the infinitely differentiable function 
q (by definition, (q ■ h,ip) = (h, q ■ if), \f(p G D as in f2R \W^), since the operation of 
multiplication is not defined for two distributions. Using this definition, it is easy to 
check that -£r(q ■ h) = h ■ -^q + q ■ -£rh. 

By proposition 3, to get v G if _1 (f2) it is enough to show (see (J27]l ) that 



qhi G L 2 (Q), 



hi-^q G L 2 (Q). 



(28) 



The first inclusion in fl28|) follows from Lemma and 



n 



\q(x)hi(x)\ 2 dx < sup{|g(x)| : x E Q} 
The second inclusion in fl28l) holds due to 



P~ 2 

fmin 



< + OO. 



a Mx) ^ q{x) 



' dx < 





9 ( \ 


sup 1 





x G tt, i — 1, no 



\hA 



< 



Pr, 



Pmin 



\h» 



< +00. 



Here we use (see Lemma) 

d 



dxi 



q{x 





9 -u \ 




dp(x) 




dx~ P {X) 




dxi P 



< 



sup 



: x G Q, i — 1, 



n 



P~ 2 

fmin 



,-2/ 



— frr, 



P~ 2 . 

fmin 



So we get the property v G H^ 1 ^) = D^A" 1 ^ 2 ) which is very important for us to 
justify the choice of assumptions on the state-dependent function £ (the choice of a 
class of functions £) in this article. Now we see that, assuming (in addition to (l25j) ) 
that f°_ r \x{9) \ dO < 00, the function £ defined by fl2T|) with v defined by fl26|) possesses 
the property ([5]). 

As a result, we may conclude that for any (finite or infinite) sequence of isolated 
points {u st ' k } C Hq(Q) C L 2 (Q) we can define a state-dependent function £, which 
satisfies assumptions of Theorem 1 and such that system ([T]) with this £ will have 
all the points {u st ' k } C Hq(£1) as stationary solutions u k (t) = u st,k ,t G R. The last 
property means that our model with distributed in space and time state-dependent 
delay term may be successfully used having information (say from experiments) on 
an arbitrary set of isolated stationary solutions. 

We notice that the definition of values of £(•, ■ , u st,k , u st ' k ) by ( J2TT) . (1261) on a set 
of isolated points does not contradict property (jl]) since the last one deals with the 
case of convergent sequence of points in H. 

To conclude this section we collect all the assumptions on functions £, b, f used in 
our considerations: 
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Ab) Function b : R — > R is locally Lipschitz, bounded and satisfies b(w) > for all 
w ^ 0. 



Afl) Function / : — Q — > i? is bounded. 



Af2) Function / is strictly positive and / G C°°(fi — fi). 

Afl) Function f : [-r, 0] x Q x L 2 (fi) x L 2 (-r, 0; L 2 (fi)) -> R satisfies the following 
condition: 

for any M > there exists £§,m so that for all (v l ,ip l ) G if satisfying ||i> l || 2 + 
S-r ll^*(*)ll 2ds <M 2 ,i = 1,2 one has 



v 1 -v 2 \\ 2 + f° \\iP\s) -ij 2 (s)\\ 2 ds 

J—r 



dQ 

1/2 



Af 2) There exists _i/ 2 ) > so that 

t r U(0,-,v^)\\D(A-^)d9 < %_ 1/a) for all G if. 

Af3) Function f satisfies £(•, ^) G L°°(—r, 0; -D(A -1 / 2 )) for all (u, ^) G if. 

Af4) There exists C (€i0 ) > so that /° r ||£(0, • ,f, ^)|| d# < <% ) for all (u, ^) G H. 

Ax) Function x satisfies /° r x(#)d# ^ and J° r |x(0)| ^ < oo. 

As an application we can consider the diffusive Nicholson's blowflies equation (see 
e.g. [261 [21]) with state-dependent delays. More precisely, we consider equation flTJ 
where —A is the Laplace operator with the Dirichlet boundary conditions, Q C R n ° 
is a bounded domain with a smooth boundary, the function / can be a constant 
as in [2S1 Ell which leads to the local in space coordinate term or, for example, 
f(s) = "^=e~ s Z 4 ", as in [25] which corresponds to the non-local term, the nonlinear 
function b is given by b(w) = p ■ we~ w . Function b is bounded and b(w) > for all 
w 7^ 0. As a result, we conclude that for any functions f satisfying conditions of 
Theorems 2 and 3 the dynamical system (St, H) has a global attractor (Theorem 3). 

So, our system (0Q) with distributed in space and time state-dependent delay term 
may be successfully used to study Nicholson's blowflies equation with an arbitrary 
set of isolated stationary solutions. 

Acknowledgements. The author wishes to thank Hans-Otto Walther for bringing 
state-dependent delay differential equations to his attention. 
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